Sudden chain energy transfer events in vibrated granular media 
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In a mixture of two species of grains of equal size but different mass, placed in a vertically vibrated shallow 
box, there is spontaneous segregation. Once the system is at least partly segregated and clusters of the heavy 
particles have formed, there are sudden peaks of the horizontal kinetic energy of the heavy particles, that is 
otherwise small. Together with the energy peaks the clusters rapidly expand and the segregation is partially 
lost. The process repeats once segregation has taken place again. Depending on the experimental or numerical 
parameters, the energy bursts can occur either randomly or with some regularity in time. An explanation for 
these events is provided based on the existence of a fixed point for an isolated particle bouncing with only 
vertical motion. The horizontal energy peaks occur when the energy stored in the vertical motion is partly 
transferred into horizontal energy through a chain reaction of collisions between heavy particles. A necessary 
condition for the observed regularity of the events is that chain reactions involve most of the heavy particles. 

PACS numbers: 45.70.-n, 45.70.Mg 



Introduction. Granular media, when externally excited, can 
have a fluid-like behavior showing many flow regimes resem- 
bling those of molecular fluids even though there are strong 
differences due to the energy dissipation at grain collisions 
111]. Fluidized granular media remain in a far from equilib- 
rium regime, presenting phenomena that makes difficult the 
development of hydrodynamic models with quantitative pre- 
dictive power. Among other specific phenomena, we remark 
absence of scale separation |2], lack of energy equipartition 
Jstl, violation of the fluctuation-dissipation relations |4], non 
thermodynamic fluxes |5], spontaneous development of inho- 
mogeneities [60, segregation |7], and localized surface struc- 
tures A careful description of these and other phenomena 
are crucial to develop models that describe the collective be- 
havior of granular fluids. 

Fluidized granular media in a shallow geometry (quasi two 
dimensional) has attracted attention because it allows for a de- 
tailed analysis of both the collective behavior and the motion 
of individual grains L9-i2J. The possibility of quantifying the 
system's dynamics at both scales may help building a mathe- 
matical model for the collective dynamics of granular media. 
Placing monodisperse inelastic spheres in a vertically vibrated 
shallow box of height less than two particles' diameters, a par- 
ticular phase separation takes place: grains form solid-like re- 
gions surrounded by fluid-like ones, having high contrasts in 
density, local order and granular temperature This phase 
separation is driven by the negative compressibility of the ef- 
fective two dimensional fluid |10]. For shallow systems the 
horizontal kinetic energy of the grains can be quite different 
from the vertical kinetic energy. 

Granular matter is usually polydisperse, with grains differ- 
ing in mass, shape, size or mechanical properties. It is known 
that a mixture of two types of grains differing in some of these 
properties can mix or segregate when externally excited \T^. 

In this letter we report an experimental and numerical study 
of a phenomenon that takes place when two particle species 
of equal size but different mass are put in a vertically vibrated 



shallow box. The focus of the present letter is not on segrega- 
tion, but on a quite peculiar phenomenon that takes place once 
the species have segregated: the horizontal energy of the sys- 
tem has sudden and brief peaks, together with fast expansions 
of the regions rich in heavy particles. This is what we call 
chain energy transfer events. The characterization of these 
events shows a direct link between the individual particle dy- 
namics and the observed collective behavior. 

In experiments, the two species do not only differ in their 
mass but, being made of different materials, they also dif- 
fer in their mechanical properties. Numerical simulations al- 
low then to identify the key elements that generate this phe- 
nomenon. Simulating the simple hard sphere model it be- 
comes evident that the mass difference is the key parame- 
ter that controls the occuiTence of energy peaks and sud- 
den expansions; the differences in mechanical properties only 
change them quantitatively. Simulations give information on 
the particle's vertical motion, helping to complete the descrip- 
tion of the phenomenon. 

System configuration. Two species of spherical grains of 
the same diameter a but different mass are placed in a square 
shallow box as described below. The lighter/heavier particles 
win be cafled L / H. 

In the experiments, we use a mixture of Nh = 170 stain- 
less steel (density 7.8g/cm^) and Nl = 291 polyamide (den- 
sity 1.14g/cm^) particles, both of a = 3 mm. The shal- 
low box has dimensions Lx/o — Ly./a = 33.33 ±0.03 and 
L,/(j = 1.813 ± 0.004. The box consists of two 10 mm thick 
ITO coated glass plates separated by a square steel frame. 
It is fixed to a base, which supports an array of high inten- 
sity LEDS. This whole setup is forced sinusoidally with an 
electromechanical vibrator. An accelerometer is fixed to the 
base of the cell, which allows the measurement of the imposed 
forcing amplitude. Top view images are obtained with a high 
speed camera. Particle positions and velocities can be deter- 
mined at sub-pixel accuracy. The vibration frequency is fixed 
to/ = a)/27r=l/r = 100 Hz. Two values of the vibration 



amplitude are used, A/cT = 0.036 and A/a = 0.045. 

In the simulations, we use instantaneous collisions that are 
characterized by restitution and friction coefficients, which 
are velocity independent and different for both kind of par- 
ticles, but the same for particle-particle and particle-wall col- 
lisions. An event driven algorithm is used |13]. The shal- 
low box oscillates vertically with a biparabolic waveform. All 
simulations start with a random initial mixing. 

Here, we report two types of simulations: The first type 
(SI) corresponds to simulations performed at higher densities 
than the experiment. Its purpose is to show that the chain en- 
ergy transfer events are robust. In these simulations, the box 
has flat bottom and top walls and horizontal periodic bound- 
ary conditions. We use the following parameters: mass ratio 
between a heavy and a light particle mn/mi = 10, box size 
Lx/ CT = Ly/a = 40, L;/ a = 1 .82, normal and tangential resti- 
tution coefficients aL~CiH = 0.8, static and dynamic friction 
coefficients iJ.s — Q.3 and /ij = 0.15, amplitude and frequency 
of vibration A/c7 ~ 0.15, co^/ajg = 7. The number of parti- 
cles are Nl — 1000 and Nh — 500, giving the filling fraction 
p = NKa^/4L,Ly = 0.736. Other box sizes, filling fractions 
and heavy particle concentration were simulated with results 
similar to the ones we are presenting. 

The second type of simulations (S2) are performed to find 
which experimental imperfections are relevant for the partic- 
ular observations (for example, the specific time and energy 
scales). Possible imperfections are the plate surface roughness 
and small curvature. The latter turns out to be more important 
than the former. In order to consider this, we use a slightly 
parabolic bottom plate, with a variation of height between the 
bottom and flat top plates of only 1%. In these simulations we 
use lateral hard walls. Comparing S2 with SI we learn which 
are the essential ingredients of the phenomenon and which 
aspects only change quantitatively the observations. For S2 
simulations, we use the following parameters: mn/mi — 7.8, 
Lx/a = Ly/a = 33, L-Ja = 1.78, ai = 0.97, an = 0.94, 
11, = 0.1 and pLci = 0.05, (O^ajg = 11.0 and A/a = 0.055 
or A/a — 0.062. The number of L and H particles are the 
same as in experiments, with p = Nna^ /AL^Ly = 0.44. 

Unless explicitly stated, the following descriptions refer to 
both simulations and experiments. Shortly after starting from 
a homogeneous configuration segregation between the two 
species takes place and many small dense clusters of //'s ap- 
pear. Later on, the clusters coalesce and tend to have some 
L's in their bulk. From a top view the //'s appear as if they 
were standing still, while the L's outside the clusters show a 
significant horizontal agitation. The external pressure exerted 
by the L's leads the //'s to form dense clusters lll4ll . The coa- 
lescence process, that is normally slow if the box is perfectly 
flat (SI), is sped up if there is a small curvature of the plates, 
coming from slight imperfections of the experimental setup, 
or imposed explicitly in a numerical simulation (S2), with the 
resulting cluster located at the center. In either case, the sub- 
sequent description of the phenomenon is independent of any 
small curvature. It only depends on the density of the cluster 
of the H's and on the number and distribution of light particles 
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FIG. 1 . (Color online) Sequence of top view configurations for simu- 
lations S 1 (top), S2 (middle) and experiments (bottom), showing the 
development of an energy peak. Heavy (H) particles are presented in 
blue and light (L) ones in red. First, the system reaches partial seg- 
regation. Then the beginning of the event is seen as a small region 
of lower density within the cluster. Next, the low density region has 
propagated to practically all the cluster. Later (not shown) the system 
comes back to a state similar to the one before the expansion. 



which remained trapped inside it. 

The chain events. Once there is at least one cluster there is a 
sudden phenomenon of short duration. Figure[T]shows typical 
temporal sequences of the grain positions when an event takes 
place, for simulations SI, S2, and experiments. A chain reac- 
tion begins as an abrupt increase in the horizontal agitation of 
the i/'s in a small region of one of the clusters, implying a lo- 
cal expansion followed by a fast propagation of the horizontal 
agitation to a larger zone. After a short while the horizon- 
tal agitation of the H's quickly decays, recovering its original 
value. Next the L-H collisions slowly compress the cluster 
again, eventually recovering the original density. Both cluster 
size and horizontal agitation, computed as mean quadratic ra- 
dius 8 and the average (per particle) horizontal energy of the 
H's, E[jh, respectively, are shown in Fig. |2]for a sequence of 
several chain events in experiments and S2 simulations. To- 
gether with the fast expansions of the cluster there are energy 
peaks. Later, 5 slowly decreases until the next energy burst 
takes place. Indeed, after the expansion, the dynamics of 5 is 
much slower than that of the energies, because the former is 
related to a conserved quantity (the mass), while energy is not 
conserved and evolves in fast temporal scales. 

Simulations show that the system is highly anisotropic, with 
vertical energies much larger than the horizontal ones for both 
species. Before an energy peak the horizontal agitation of the 
L's is higher than for the H's. The reason is that the L's are 
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FIG. 2. Time evolution of Em, (a), S (b) and alf (c) for an experi- 
ment witfi A/a = 0.045. (d), S (e) and cThz (f) for a simulation 
SI andA/cT = 0.062. 



having collisions with particles of both types and, in particu- 
lar, the collisions with the H's keep them excited both verti- 
cally and horizontally. On the other hand most i/'s are at their 
fixed point, with almost all their energy concentrated in the 
vertical motion. The abrupt change of £///, is accompanied 
by simultaneous but weaker changes in the other computed 
energies. This behavior is generic. 

Under certain conditions, a single grain bouncing in a shal- 
low box vibrating at a frequency CO evolves towards a fixed 
point: the bouncing movement becomes periodic with the 
same frequency co and, because of the friction with the hor- 
izontal walls, the particle does not move horizontally neither 
does it rotate. Namely all its kinetic energy is vertical energy. 
The role of the fixed point in the behavior of the whole system 
is evident by looking at side view videos of experiments and 
simulations 1141 . It is seen that the H's are bouncing in phase 
as practically one solid layer: they are collectively trapped at 
the fixed point very close to each other. In fact, the fixed point 
can be easily computed in absence of gravity (a good approx- 
imation in this case) with a resulting vertical energy for the 
heavy particles in the simulations of 0.32 mi(aa))^, that fully 
agrees with the value of Ehv before and after the energy-peak 
event, confirming that most H's are close to the fixed point. 
Hence, the collisions among heavy particles create almost no 
horizontal agitation. It is this coherent movement that is de- 
stroyed when an energy-peak event takes place. 

To check that the picture described above is correct we 
have measured in simulations the standard deviation Ohz of 



the stroboscopic height of the H's when the box is at its low- 
est position, shown in Fig. |2f). Ohz is almost zero before 
the energy burst and it jumps an order of magnitude during 
the event. It recovers its typical small value together with 
Ef{h^ indicating that the H's are again moving in phase. The 
cross correlation function between Efji, and Ghz shows a clear 
maximum centered at null time delay, confirming that the en- 
ergy bursts are accompanied by a massive desynchronization 
of the H's. Experimentally, this sudden desynchronization 
is observed through aif, which is the envelope of the ac- 
celeration signal including particle-plate collisions, obtained 
by low-pass frequency filtering. The main contribution to 
flLF comes from the simultaneous H particle-plate collisions, 
which therefore exert a large force on the plate. The negative 
peaks, that occur at the same instants than the energy peaks 
(see Fig. |2|c)), correspond to the absence of coherent colli- 
sions due to the dephasing of the H's and the decrease in their 
vertical kinetic energy. 

In summary, at the start of the dynamics the horizontal 
energy £///, of the H's begins to decrease, reaching a small 
asymptotic value because most of the H's are approaching the 
fixed point. Particles L produce only small perturbations to 
the heavier ones but, due to the mass difference, the H's can 
easily take the L's off their fixed point getting significant hor- 
izontal energy. As a result, the L's very seldom reach their 
fixed point and if they do so, it is for a short time. The per- 
manent collisions with the L's that some of the H's suffer, 
allow for a residual non vanishing horizontal energy Ehi,- The 
clusters of H's slowly evolve to higher densities. At one mo- 
ment, an H may collide with an L in such a way that the H 
gets a large horizontal velocity. The heavy particle then hits 
a neighboring energetic H triggering a chain reaction of colli- 
sions among neighboring H's rapidly propagating within the 
cluster implying a sharp peak in E[fh ■ The main energy source 
for the increase of E^h is the vertical energy Ehv which cor- 
respondingly decreases during the energy burst. The chain 
energy transfer induces that the En, and L^,, follow a similar 
pattern to a smaller extent. This sudden chain energy transfer 
from the vertical to the horizontal motion is responsible for 
the observed fast increase of the cluster radius 5. 

Statistical analysis. The time lapses between successive 
Efji, peaks and their intensities show different degrees of reg- 
ularity; in some cases the peaks seem almost periodic. Fig. [3] 
presents the power spectral densities of 5/(7 for both experi- 
ments and simulations S2. Two distinct regimes are observed: 
for low amplitude there is no periodicity in the signal, while 
for A/(7 = 0.045 for experiments (A/ a = 0.062 for S2 simu- 
lations) a clear maximum, although broad, is observed. It is 
centered at a period « lOOOL (« 400r for simulations), which 
is consistent with the time series in Fig. ^ (Fig. |2^). Depend- 
ing on the system and forcing parameters, simulations S 1 also 
show almost regular regimes. In some cases these are much 
more regular than those obtained in S2 simulations. 

After exploring many different system parameters there ap- 
pears to be an explanation for the observed periodicity un- 
der certain conditions. We observe that in systems where the 



4 




0.1 0.2 0.3 0.4 0.1 0,2 0,3 0,4 

f (Hz) f (Hz) 



FIG. 3. (Color online) Power spectrum density of d/o. (a) Ex- 
perimental results for A/cT = 0.036 and A/<7 = 0.045 (blue and red 
curves respectively); (b) Simulational results (S2) for A/a = 0.055 
and A/<T = 0.062 (blue and red curves respectively). 

region of horizontally excited H\ always propagate through 
most of a unique cluster, the peaks are more regular This is 
so because for the chain reaction to take place, a high enough 
density of the cluster is needed, and it takes a characteristic 
time for the cluster to reach this density after such an event. 
If the chain reaction does not propagate through all the cluster 
then the time of compression is highly variable as it depends 
on the amount of particles that were involved in it, and also 
other similar events can take place while the cluster is com- 
pressing in those parts of the cluster that remained dense. In 
order for the explosion to cover most of the H cluster, sev- 
eral conditions must be fulfilled, particularly that its density 
is large enough to allow for an uninterrupted propagation of 
the horizontal energy. Moreover, the concentration of L's in- 
side the cluster needs to be small so that they do not block the 
propagating front. These conditions can be achieved either by 
changing the control parameters or by having small geometri- 
cal defects such as a small curvature of the box to help the //'s 
migrate toward the center Indeed, in S2 simulations, an arti- 
ficial small curvature was added and the resulting events are 
more regular compared to the same simulations done without 
the small curvature. 

Conclusions. Experiments and numerical simulations of a 
granular mixture of grains that differ in their mass density — in 
a vibrated shallow box — show horizontal energy peaks char- 
acterized by the rapid conversion of vertical kinetic energy 
into horizontal one; these events are preceded by the segrega- 
tion of the species. In the segregated state the massive grains 
approach a fixed point characterized by a vanishing horizon- 
tal energy and vertical motion in phase with the walls' oscil- 
lation. The synchronized heavy grains collide between them- 
selves with a tiny momentum transfer. As a result, densely 
packed clusters of the heavy grains develop. Eventually, how- 
ever, a light grain can interpose in the vertical motion of a 
heavy grain and, consequently, the latter acquires significant 
horizontal momentum and leaves the fixed point. This can 
also be triggered by another type of strong fluctuation, like 
a sudden desynchronization of one H particle due to surface 
roughness or an energetic incoming L particle. The subse- 
quent collisions with neighboring heavy grains transfer a large 



amount of energy into the horizontal agitation in the form of 
a chain reaction, generating an £///, peak and an expansion 
of the heavy particles cluster When there is a unique clus- 
ter different regimes are reached depending on the density of 
the cluster of heavy particles and the relative concentration of 
light particles inside it. When there are too many light par- 
ticles in the cluster the horizontal agitation does not involve 
the entire cluster and the corresponding chain events show no 
characteristic time or energy scales. Otherwise, explosions 
propagate through all the cluster and show a characteristic in- 
tensities and time lapse between successive events. 

Simulations and experiments differ in the collision details 
and, consequently they do not give the same quantitative re- 
sults. Although the fixed point exists only for spherical par- 
ticles, small non-sphericity, as present in the experiments and 
artificially put in simulations, result in a small horizontal agi- 
tation but the energy bursts still take place. The independence 
of the phenomenon on the precise details shows that it is ro- 
bust, and requires only a large enough mass contrast. This 
phenomenon shows that to correctly describe the collective 
dynamics in a confined geometry, the dynamics of individual 
grains in the vertical direction is crucial. The revealed link be- 
tween the small and large scale dynamics should help building 
hydrodynamic-like models of confined granular media. 
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